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ABSTIUCT 



Tlie objGctive of the thesis Is to minimize the total 
time required for a driver to pass through the Lincoln 
Tunnel, travel to his destination in the Central Business 
District of New York City, wait, be serviced and then exit 
his destination. This is accomplished by requiring vehicles 
to delay in New Jersey before entering the tunnel, 

A general model for total time in the system as a 
^ function of the waiting time in New Jersey is presented. A 
sampling teclinique is developed vdiich tests the adequacy of 
that portion of model which predicts the conditions within 
New York City. 

The control techniques are extended to apply to 
ships entering harbor through the constrains of a channel. 
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CHAPTER 1 



DESCRIPTION OF THE SYSTEM 



I. 1. Introduction 

The central business district (C.B.D.) of New York 
City, south from 42nd Street, may be viewed as a large and 
complex commerical and industrial center. Over the years, 
activities of a similar nature have tended to locate 
together and thus somewhat separate, although not entirely 
independent, business districts have formed. 

There are numerous entrances into and exits out of 
the C.B.D. From the point of view of the vehicles travel- 
ing into New York City, the various business districts may 
be thought of as traffic receptors or traffic sinks. Hence- 
forth in this thesis, these will be referred to as traffic 
sinks or simply as sinks. At the present time there is no 
attempt made to limit the number of vehicles that enter the 
C.B.D. and, as a result, large traffic congestions are 
common. Consequently, a driver entering New York City can 
expect long delays before concluding his business. 

The attempt of this thesis is to reduce these delays 
by requiring vehicles to wait before entering the C.B.D. 

This will allow the particular traffic sink to clear out 
enough so that the total time a driver spends in New York 
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City will be reduced. One entrance, the Lincoln Tunnel, 
will be studied in detail. 

1.2. Business Districts in New Yoric City 

Although no rigorous attempt has been made to 
delineate the New York City traffic sinks, the Port of New 
York Authority has supplied the author with sufficient infor- 
mation to describe roughly the boundaries as shown in Fig. 1. 

The author attempted to determine exactly the limits 
of one of these sinks. It was found that vdiile difficult 
and time consuming, such an undertaking was not impossible 
and could be effected if sufficient personnel and time vrere 
available. Fig. I is time dependent; for example, districts 
12 and 13 are in the process of relocating further north. 
Other districts can be expected to expand, contract or 
relocate as time goes on; this indicates continuous updat- 
ing. 

1.3. The Lincoln Tunnel 

Some of the entrances to the New Yorfc City C.B.D. 
are also shown in Fig. 1. The Lincoln Tunnel is a three 
tube, six lane facility operated by the Port of New York 
Authority and represents a cost of approximately $18,000,000 
per lane per mile.^^^ Each tube has 2 lanes and is about 
7,650 feet long. The south tube accommodates east bound 
traffic exclusively, the north tube handles only west bound 
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Government Complex 
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Financial District 
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Import/ export District 
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Fish Market 


5 


Electronics (Wholesale and Retail) 
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Jewelry District 
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Residential 
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Textile Manufacturing 
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Offices 
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Light Manufacturing 
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Greenwich Village 
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Wholesale Meat Market 
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Produce Market 


14 


Residential 


15 


Mens Clothing 


16 


Garment District 


17 


Residential 


18 


Hospital Center 


19 


Residential 


20 


Insurance Center 


21 


Toy District 


22 


Air Line District 


23 


Offices 


24 


Residential 




FIGURE 1 

The Traffic Sinks of New York City 
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traffic, and the center tube is switchablo. During the 
morning rush hour traffic, the center tube is used to trans- 
port vehicles into New York and during the evening rush hour 
it is changed into a west bound roadway. At other times the 
central tube handles traffic flows in both directions. 

Although rather complete information on the origins 
and destinations of vehicles vAxich use the Lincoln Txannel 
has been collected^^^^ and data on the behavior of vehicles 
within the tunnel has been obtained, (26) j^q success has been 
realized in attempting to optimize the use of the center 
tube. 

Since the adequacy of the control techniques, to 
be presented in later chapters, depends upon the traffic 
engineer's ability to control the traffic within the 
Lincoln Tunnel, some of the theories of traffic flow will 
be examined in the next section. 

1.4 General Theories of Traffic Flow 

Many authors have studied traffic flow in detail and 
only a brief review will be presented here. One relationship 
basic to the study of traffic flow is response at time (t+T) 
equals sensitivity multiplied by stimulus at time (t).(^^^ 
Response is usually taken as acceleration since that is the 
variable under the driver's control. The stimulus is gener- 
ally represented as the difference in the velocities of two 
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vehicles. Thus, ^ (X^(t+T))* sensitivity x minus 

dt 



dt n+J- 



(t)) where X represents distance traveled. For 



mathematical simplicity, sensitivity was originally taken as 
a constant. Later, however, more sophisticated models of 
sensitivity have been developed^ 

A fundamental diagram of traffic flow has been 
developed^ and is shown in Fig. 2, It is stated that 
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FIGURE 2. 

Fundamental Diagram of Traffic Flow 
such a diagram is valid only for a particular section of 
roadway at a particular time. Authors have placed different 
interpretations on this fundamental diagram, Martin^^^^ 
feels that it is valid only in equilibrium or steady state 
conditions and as such represents the desires of the drivers, 
Duckstein believes that Fig. 2 is the center of a band of 
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possible flow-density combinations. Many experimenters feel 
that traffic density (k) is the best indication of the 
actual traffic conditions^ Average velocity is not a 
good indication because small speed changes can represent 
large changes in flow, Greenberg has compared traffic 
flow to that of a continuous fluid, 

Shock waves appear to be a major cause of low flow 
rates, It has been shown that a small increase in peak 
flow rates results in dramatic reductions in traffic conges- 
tion. <^0) 

Authors feel that density can be used to predict 
the formation of a shock wave^^^ and that high output speeds 
indicate the absence of a shock wave. Officials of the 
Tunnels and Bridges Department of the Port of New York 
Authority have demonstrated that there exists a critical 
output velocity and that traffic densities can be used to 
foresee sub-critical velocities. ^ Finally, it has been 
shown that inherent Instability exists in non-homogeneous 
traffic flows, The velocity at \^ich traffic actually 
travels is a compromise between the desired velocities of 
various components of a non-homogeneous flow. Thus the 
components which are required to travel below their optimum 
velocity can react to a sudden decrease In the speed of 
leading vehicles In such a way as to create a shock wave. 
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The Port of New Yoric Authority has been and is 
presently attempting to induce small, frequent gaps in the 
traffic flow. These small spaces do not decay as rapidly 
as larger, less frequently spaced gaps. Only a small degree 
of success has been realized in implenenting these goals. 

The south tube of the Lincoln Tunnel is equipped with a 
closed circuit television monitoring system, but as of this 
writing, the tunnel officials only react to but do not pre- 
vent the formation of shock waves. 

The traffic sinks in New York City and the Lincoln 
Tunnel have been discussed briefly. This is necessary for 
an xinderstanding of the objective of the thesis which is 
presented in the next section, 

1.5, Objectives of the Thesis 

Once a vehicle departs the ttinnel, his problems have 
just begun, A driver will find himself fighting traffic 
all the way to his particular traffic sink and then fighting 
for a position in the sink. It is felt that by requiring 
a vehicle to wait in New Jersey before entering the tvinnel 
it may be possible to allow his traffic sink to depopulate 
sufficiently so that the total time required for him to be 
serviced in New York, is actually reduced. The basic idea is 
to ranove this driver from the congested area long enough 
for some vehicles in the area to depart. 
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Since it is almost an impossible task to determine 
the destinations of cars and small vehicles, it was decided 
to apply this principle to large trucks and tractor trailers 
which supply and remove goods from the C.B.D. of New York 
City, 

The general model is presented in Chapter 2, Two 
basic probl«n areas arise; (1) sampling the sinks to 
estimate various parameters and (2) determining if the 
tunnel supplies a large enough percentage of the sink 
traffic to make effective control of the sink conditions 
possible. These problems are treated separately in Chapters 
3 and 4, 



CHAPTER 2 



MODELING THE SYSTEM 



2.1, General Model 

As a vehicle approaches the Lincoln Tunnel In New 
Jersey, It enters a large and complex traffic system. This 
system consists of the waiting areas in New Jersey, the 
tunnel, the network of roads in New York leading to the 
destination, and the traffic sink itself* Upon being ser- 
viced and departing from the traffic sink, a vehicle is con- 
sidered to have exited the system. The Lincoln Txannel may 
be thought of as an entrance way leading to 24 separate but 
not wholly independent sub-systems or traffic sinks. As a 
vehicle leaves its sink it may travel into another sink 
and thereby become part of the xmcontro liable input into 
another sub-system 

Total time in the system (T) consists of six compon- 
ents (t.) such that ,E,t. =» T where: 

t^ s waiting time in New Jersey, 

t 2 = travel time in the txannel, 

t^ * travel time from tunnel exit to desired traffic 
sink, 

t^ = waiting time in sink. 
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= service time in sink, 

■ time required to exit sink after being serviced. 
The time required for a vehicle to return to its 
original destination is assumed constant and is not con- 
sidered in the model. All times are measured in minutes. 

This model is valid for any time period during the day, how- 
ever, here it will be applied to the morning rush hour 

traffic between the tLmes t ^ ^ and t See Fig. A. 

start end 

All times are assximed to be between these limits 

2,2. Components of the Model 

The individual components of T will now be investi- 
gated separately. 

2.2.1, Waiting Time in New Jersey 

Waiting time in New Jersey (t^) is the independent 
Variable. T will be minimized by controlling tj^. An upper 
bound of 25 minutes has been chosen for this variable, 

2.2.2, Travel Time in the Tunnel 

It was felt originally that in order to predict t^ 
satisfactorily it would be necessary for the model to predict 
the status of each vehicle in the tunnel at every second. 

This involved estimating distance traveled, instantaneous 
velocity and acceleration and headway for every v^icle. It 
soon became obvious that this approach would be of little 
practical value because of the large amount of computer 
memory space and computing time required to predict t 2 » 
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To avoid these problems, it was decided to utilize 
the f low-d'ensity data available for the tunnel, A quad- 

ratic curve was fitted to this data by the method of least 
squares, and the solution of the normal eqxiations was; 

0 a a + bk + ck^ where 
0 = flow measured in vehicles/minute, 
k = concentration measured in vehtcles/mile, 
a = 0,44, b = 0.61, and c =-0,0043, 

The predicting curve is shown in Fig, 3, 




FIGURE 3 

Flow (e) vs Concentration (k) 
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It Is also possible to observe flow as a function of 
time of day. This data was obtained and the flow between 
the times of tg^art ^end approximated by a quadratic 

function. This curve is shown in Fig. 4, The quadratic 




FIGURE 4. 



Flow (e) vs Time of Day (t) 

equation is 0 * d t^ + pt q with d = -9.2, p = 127.8 and 
q = -241. 

If it is assiimed that no shock waves will appear in 
the traffic flow, then 9 remains fairly constant for any 
given time of transit through the tunnel. For any time of 
the day, 

2 

ck '►bk + a-esO, 

-b^^/b^«»4c( a-9 ) 



k 



2c 
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1' = ~-6lj:^.61^»A(..0043)(.44»9) 

2(-.0043) 

For t < t < t .,6 will always be greater than 

start - - end* 

.44 but the expression under the radical will always be 
positive. This can be seen by substituting for G, 

Since the model predicts traffic concentration for rush 
hour traffic, it seems reasonable to asstime that the actual 
density is the larger of the two densities corresponding to 
a specific value of G, Thus, 

k . 6 iW.6i'^-4(-.0043’) ( . 44-9^ 

2(.0043). 

The average velocity through the tunnel (v) is the flow 
divided by the density. Thus, velocity in miles/minute 
may be written as; 

V == 9/k, 

V = 29c/(-bV^b^-4c(a-9)). (1) 

As expected, ^^q. 1 simply says that on the heavy concentra- 
tion half of Fig, 4, velocity increases then decreases with 

e. 

o 

It was shown that 9 = dt + pt + q for ^ £ 

^end expression is substituted into Eq. 1 to obtain 

velocity as a function of time; 

V =s 2c(dt^+pt+q)/ (-bWb^-4c(a-dt^-pt-q) ) , 

Since this velocity is considered to be the average velocity 
through the tunnel, the expected time required for a passage 
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through tho tunnel is found by dividing the tunnel length 
(L miles) by this average velocity. Thus, after inserting 
the values of the parameters, 

= Lt ■ 0043) (~.iW.T6t^-33t*97 ) ) ( 2 ) 

(-2. 36t^+33t.97) (2) (.0043) 

Since t, the time at which a vehicle enters the tvmnel, is 
the time of arrival at the New Jersey entrance plus the 
waiting time In New Jersey, it is seen that t 2 depends on 
the control variable tj^. As t varies between and 

^end* ^ increases and then decreases (Fig. 4). This results 
in the unlikely conclusions that t 2 reaches its minimum 
Value when t = t*, or when 0 is at its maximum value. This 
is shown graphically in Fig. 5 in which t 2 is plotted as 
a function of t as t varies between and 

This erroneous result was cuased by the two follow- 
ing unrealistic assvimptions: 

1, The actual concentration is always tho larger 
of the two concentrations corresponding to a specific 9, 

2, No shock waves will be induced in the flow. 
Although It is true that during peak traffic the flow is on 
the heavy concentration side of Fig. 2, it is not reasonable 
to assume that the flow started on that side. As the morning 
rush hour approaches, flow increases along the left hand side 
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FIGURE 5. 

Time in Tunnel vs Time of Day (t) 

of Fig. 3. At some point, the flow shifts over to the right 
hand side of the curve; this can happen in one of two 
ways ; 

1, Flow can reach its maximum from the left and 
then decrease down the right side of Fig, 3, 

2. As flow is approaching its maximum value, a 
shock wave can shift the flow to the heavy concentration 
side of the graph. Thus the true expression for t 2 as a 
function of t is; 
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t = L(. 0043k. 4&t9.2t'‘ ^. 127. 8tt42l)) 

2 (-9.2t“+127.8t-421)(2)(.0043) start- I 

or 

t = i-l. 61 1- ' . 6 1^.4( . . 0043) 2t^- 127 . 8 tt42 1) ) 

" (-9.2t^+L27.8t-42L)(2)(.0043) 

I 

whore t represents the time at which maximum flow is 
achieved if no shock wave occurs or the time of a shock wave. 

The solid lane in Fig. 6 shows how t 2 varies with t 
if no shcok vrave is induced; the dotted line traces t 2 if a 

t 

shoclc wave is encountered at t. In either case, it can be 
concluded the t 2 is a non-decreasing fxinction of t. 

If Fig, 4 were used as a guide, one would decide 
that the flow will never bo on the light concentration side 
of the flow-density curve until evening rush hour. This 
obviously is not necessarily true. An acceleration shock 
wave can occur as 9 is decreasing along the right side of 
Fig. 3. If this happens, t 2 will begin to decrease. Fig. 

7 shows this penomena graphically. 

2.2.3 Travel Time to the Traffic Sink 
The author had hoped to study t^ in detail, but 
because of time limitations, this was impossible. Therefore 
t^ will be taken as a constant for each traffic sink. Refer- 
ences (2) and (23) among others, discuss some of the prob- 
lems that would be encountered in predicting t^. 



Time in Tunnel 

(i-iijiutes) Time in Ttinnel (Minute 
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FIGUP.n 6 . 

Time in Tunnel (12) vs Time of Day (t) 







FIGURE 7. 

Time in Tunnel (12) vs Time of Day (t) 
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2.2.4 Time in the Traffic Sink 

It is almost impossible to know exactly vhat will 
happen to a vehicle once it enters a traffic sink. Some 
vehicles will find a legitimate servicing station immediate- 
ly, some will create their own servicing stations by double 
parking or by parking in a restricted zone, and some 
vehicles will be forced to wait until a servicing station 
becomes available. It is also impossible to determine 
exactly what will happen to a vehicle once it locates and 
enters a servicing station. Some vehicles will find men 
available to unload or load cargo v*iile other vehicles will 
be required to wait for servicers. As the congestion mounts, 
servicers will begin to interfere with each other and will 
increase service time. It becomes clear that the standard 
queueing models are not adequate to describe the situation. 

It is hypothesized that the total time spent in a 
traffic sink (t^ + t^ *»■ t^) increases with the number of 
vehicles in the sink at the time of arrival (N), It is 
further hjrpothesized that this increase is not linear but 
that the rate of change of the total time in a sinlc with 
respect to the number in the sink also Increases with K. 

Letting w “ t^ + t^ + t^, the following model will 
be hypothesized: E(w) » ae^^\ where a and 3 are parameters 

vrfiich must bo estimated for each traffic sink. The possi- 
bility of using other models for w is discussed in Chapter 3. 
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For the rnoment, it is assumed that t 2 = = 0 so 

that Vehicles enter their sinks directly and instantaneously 
from the New Jersey waiting areas. Total time in the sink 
will now be studied as a function of tj^. 

If it is assumed that all entries into a sink occur 
instantaneously at the beginning of each minute, then a 
consequence of the model is that all the vehicles that enter 
a particular traffic sink during the same minute have the 
same expected tLme in the sink F(w), Trucks vrhich enter a 
traffic sink during the same minute are said to comprise a 
block of vehicles. At time t, a givcsn sink may be represent- 
ed as in Fig. 8. Here K. represents the number of vehicles 
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FIGURE 8 

Conditions within a Traffic 
Sink. 



expected to depart at time t.. The expected number of 

entries per minute from all sources other than the tunnel 

is denoted by A.. It is assumed that \ remains constant 

for and that the first in first out rule is 

k 

valid for the traffic sinks. Letting N = E Ni the 

i=sa ^ 

following relationships can be devloped: 
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If = 0, E(w) = 

= 1. E(w) = 

and in general E(w) = ae^ for 

The sum tj^ + E(w) = (tj^-fl)!**-!) j.gg^,j^Qg minimum 

Value when tj^ = 0, 

For t - t<t. <t. - t The general expressions for 

B "*1 D 

expected time in the sink is E(w) = oe^ The 
sum of t|^ and w is minimized when ** Similarly 

for t^ - t < tj^ < tj - t the total time in the system is 
tj^+ae^ and is minimized when t^ = 

- t. 

It is also seen that the value of tj^ which minimizes 
T is such that a vehicle will arrive at the traffic sLnk at 
the expected time of departure of a block of vehicles and 
never between the expected times of departues of two blocks. 
Unfortunately t 2 ^ 0 and t^ 0 and it is not poss- 
ible to control tj^ so that a vehicle will enter a traffic sink 
just as a block of vehicles is departing. Considering that 
vehicles must go through the tunnel before entering their 
desired sinks, it is clear that the present expression for 
w is not adequate. When a vehicle arrives at the tunnel, the 
control personnel roust predict the conditions at his desired 
traffic sink for the coming 25 minutes. Then a decision can 
be made as to how long to delay the vehicle in New Jersey, 
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This prediction must consider the vehicles bound for this 

specific sink which have already entered the tunnel but have 

not yet arrived at the sink. 

Let M be the number of these Vehicles which have 
o 

passed the control personnel but have not yet arrived at 

the sink. Then the expression for E(w) must include a teirm 

f(M ,t,t,) which predicts the arrival of these vehicles. 
o I 

Let g(N,t,tj^) predict the dissipation of vehicles from the 
sink. The expression for E(w) is now 

E(w) = aexp(e(M^-g(N,t,tj^)+f(M^,t,tj^)+(tj^+l)X-»-l)). 



2.3, Minimizing, Total Time in the System 

Combining all parts of the model yields: T( t , t j^)=t 

L ( . 61-r^. 61^-4(.7 o 04^ ( '. 44+9. 2( t+t j^) ‘^^127'. 8 ( t+T ) •*•421) ) 

— -2 

(-9.2(t+tj^) •fl27.8(t+ti^).42l(2)(.0043) 



+ constant + 

0(N-ga;.t.tj^)-ff(M^,t,tj^)+(tj^+l)X+l) 



ae 



(3) 



where. 



t = time vehicle arrives at the tunnel 

(-1, { < (t+t^) < 

Chapter 3 discusses ways of estimating a and 3, 
Assuming for now that this has been done, there are at least 
two ways of finding the value of tj^ which minimizes T. 
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1. By assioming the form of f (M^t, tj^)and g(N,t,tj^); 

it is then possible to solve Eq. 3 numerically for the opti- 
mum value of An example of this is done in Chapter 4, 

2, By computer simulation; the case for this method 
is strong. The tunnel-sink system is extremely 
complex and it is impossible to describe mathematically 
exactly what will happen to a vehicle which enters the 
system. With simulation it is not necessary to guess f and 
g, nor is It necessary to assume first in, first out. As a 
Vehicle enters the tunnel, its expected time in the tiinnel 
is calculated. When this predicted time is reached, the 
vehicle is simply transferred out of the tunnel into the net- 
work of roads leading to the sink and its time of. arrival at 
the sink is predicted. It has been assumed that a decision 
is made for only one vehicle at a time. With simulation it 
is possible to determine the effect of letting any number of 
Vehicles through to the same sink. Time of arrival at the 
tunnel (t) is really a random variable. With simulation, 
the modeler can study the effect of variations in t as well 
as variations in the initial conditions in the tunnel and in 
the sinks. Such a program has been written. The program 
listing is in Appendix A and the logic flow chart is in 
Appendix B, 
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Before the experimenter can apply the theories of 
this chapter to a particular system, he must bo satisfied 
that he can predict the expected waiting time in the traffic 
sinks. This problaii is discussed in the next chapter. 



CHAPTER 3 



DRAWING SAI'IPLES FROM THE TRAFFIC SINKS 
3.1, .^Exponential Hodel 

As vas discussed in Chapter 2, it is assumed that 
the sum of waiting time in a traffic sink, som^icing time 
and exit ti.me (w) is equal to ae^*\ For any specific 
traffic sink, the experimenter must be able to estimate a 
and g and obtain some idea of how well his model represents 
the conditions of the physical sink. The purpose of this 
chapter, then, is to present a sampling technique by which 
these aims can be achieved. 

To be specific, the model for total time in the sink 
is : 

Wj^ = ae ■*’ where 6 i = 1, 2, ••• m, 

and m = the number of observations, 
it is assumed that: 

a, deviations are normally distributed, 

b, deviations have a common variance 6^, and 

c, deviations are independent. 

The actual graph of w vs N depends, of course, upon 
the specific values of a and p, but in general, it can be 
represented as in Fig, 9, 
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Humber of Vehicles In Traffic Sink 



FIGURE 9, 

Time in Traffic Sink (w) vs Number of Vehicles 
in Traffic Sink (N) 



Since N(0,6^) = . 1 . . exp(-V(w-ae^^)^*6"^) , the Like- 

J7rrfi 

lihood function is L=(-^)rr./2.-^* exp(- - .S. (w i) ".5*’“^) . 

2n 6 1=1 i 



Since II and are observable values, there 
meters which must be estimr tad : a and g, and 



are three para- 
2 

6 , L is maxi- 



mized by selecting a and g such that: 




w e 
i 



gKi 



a 




*2eKi = 0 



(4) 



and 



E w.N.e -^-a 
ta.1 i i 



rn 

E 

i=l 



N 



2^Nt , 



= 0 



(5) 
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Solvin5 Eqs. 4 and 5 simultaneously yields the 
maximum likelihood estimates of a and 0. It is interesting 
to note that these are identical to the lea st«> squares 
estimates of a and 0. The sum of the squared deviations of 
the observed values of from the predicted values is 
S(a,0) = (w^-ae^ , Upon taking the partial derivi- 

tivies of S with respect of a and 0 and equating them to 
zero to minimize S, Eqs, 4 and 5 again result. 

To obtain the maximum likelihood estimate of 6^ one 

9 

simply maximizes the partial of L with respect of 6 . The 

® 2 ^ R ^ ^ 

result is 6 « Z (v/. where m is the number of 

issl *• 

(wj^,Ni) tvx>- tuples observed and • over a parameter represents 
an estimate of that parameter. It is not possible to deter- 
mine whether or not ^ is biased; this requires a knowledge 
of the bias in a and Kovrever, maximum likelihood 
estimators are, in general, not unbiased. 

One advantage of maximum likelihood estimators is 
that, tinder fairly general conditions, it can be shown that 
they are efficient estimators. From this it follows 
that maximum likelikhood estimators are also consistent. 

3. 2. Determination of Sample Size 

At this point it would be desirable to have a back- 
log of samples from which it would be possible to obtain an 
advance estimate of experimental error (6^). From this the 
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experimerjtor could determine the sample size required to 
satisfy his precision raqulrments. If such a badclog of 
samples is not available; it will be necessary to obtain an 
advance estimate of 6^ in another way. 

The proposed methods of drawing samples and estimat- 
ing parameters will be demonstrated by means of an example. 
The traffic sink selected was a Circle K food market in 
Tucson, Arizona, It is not hypothesized that the Circle i' 
represents a New York traffic sink; however, the food mar- 
ket does possess some quantities similar to those of a 
traffic sink. First, most people enter Circle K with 
similar missions — to buy one or two items. In New Yori<, 
trucks entering traffic sinks have similar missions also. 
Second, there is reason to believe that the congestion 
created by the presence of many cars in a Circle K park- 
ing lot will cause the rate of change of waiting time to 
increase with N, 

For every entering vehicle, it is necessary to 
determine Wj^ and A convenience sample of size m«8 was 

drawn and the results are shown in TabLo 1, 

A modified steepest descent method using an itera- 
tive search was used to solve the nonlinear equations 1 and 
2. The following results were obtained: 

a = 0.726, p =! 0.64 and 6^ = 0,190, 
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TABLE 1. Results of Pilot Sample 



i 




(minutes) 




1 


0 


1 




2 


1 


2 




3 


2 


3, 2, 2, 3 




4 


3 


5 . 5 





How, using this estimate of experimental error, it 
is possible to determine the sample size necessary to satis- 
fy the precision requirements. These requirements may be 
worded In several ways, A few examples will be given. 

3.2,1. Requir«nent I: 

Determine the sample size (m) , such that E(w) can be 
estimated within \ minute 99,5% of the time. 

I 

Since i s d =» % minute, 



(t.995.m-l)2(i2) 

m = = (t 



. 995, m-1 



)^(0.76). 



It is necessary to iterate on m until the equality 
is satisfied. It was found that with m=8, (t ^)^(.076) 
equals 9.75, With m=9, (t gg^ 0)(O.76) = 8.51. Therefore, 
for these precision requirements a sample size of 9 is 



necessary. 



3,2.2 P.equ ir eir.en t II: 

Determine the sample size such that C ts 99% confi- 
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dence limits on E(w) will have, with probability p=.99, a 
length less than or equal to one minute. 

Here we use the relationship developed by Harris, 
Harvitz, and Hood:^^'^^ 



d2 






where 



v+1 = pilot sample size = 7, and 
d » I minute. 



Hence, 1.32 = .7^ . 

m 

Again it is necessary to Iterate an m until the equality is 
satisfied. In this example, a sample size of 35 will fulfill 
the precision requirements. 

3.2,3, Requirement III; 

Determine the sample size required for the probabil- 
ity of a 90% confidence interval on E(w) having a length 
d<6/5 to be equal to 95%. 

Letting be the estimate of £(wj^) it follows that; 
P( w£-£.'(wj^)) < = 0.90 and 

P( Wj^-E(wj^)) < 6/10) = 0.95. 



St 

— ^ < 6 / 10 . 

ym “ 



Hence, 
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lo find the mintimain value of m that satisfies the given 
requirements, set = 6/10, 

Squaring both sides and simplying, Eq, 6 is developed: 
:2 m2 



ns 



lOOt^ 



0.95,m-l* 



(6) 



? 9 

It is now observed that 



m-l 



For Values of v greater than 30, it is possible to 
use the fact that - V2V-l) 2, is distributed 

approximately as N(C,l) to obtain the following approximation, 

y-l° 

Now, incorporating this result into Zq. 6 the result is, 
'%%■! V ' ^^^.93 rni_^ 



X 



,95, m-l 



^°^^‘,95,m-l ’ 

= %(z_g5+ys^)2 = s2-^ 



loor 



,95, m-l ’ 



hence 



m 2= U + ,/2ni.3)2. 

For values of m greater than 30, the approximation 



'.95,n,-l~^-.95 

Therefore, = 50Z",g5(Z •V2n\-3)^ and because a negative 

sample size is meaningless, 

m =5 V^50 2 (2 + v*2m- 3 ) , 

m^ - 2(19, I5)m + (19. 15)^ = (ll, 62) 2(2m-3) , 

- 310, 3m + 775 +■ 0, 
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Therefore: 

3 WaTo. 32.4(1) (yTsJ 
2 

^ - 310.3 ± 303 

2 . 

n = 305 or 4. Upon substituting these two sample sizes 
into the precision requirements it is found that rn = 4 is 
too small and m = 305 is the only acceptable solution. 

It is readily seen that Requirement III is by far 
the most restrictive; it does not require an advance esti- 
mate of 5 , It is felt, however, that few people would word 
their precision requirements in this manner. Therefore, for 
this example. Requirement II was selected and a systematic 
sample with a random start of size m = 35 was drawn. Ths 
results are summarized in Table 2, 



TABLi: 2, Results of Sample 



I 


"l 


f^’i 




”i 

.Z, w, 
J=l 1 


1 


8 


0 


1.5, 1.0, 0.5, 1.0, 2.5, 1.0, 1.0, 1.5 


10.0 


2 


8 


1 


2.0, 2.0, 1.5, 3.0, 1.5, 2.0, 1.0, 2.0 


15.0 


3 


9 


2 


2.0, 2.5, 2.5, 2.0, 2.0, 3.0, 2.0, 2.0, 3.0 


21.0 


4 


7 


3 


3.5, 3.5, 4.0, 4.5, 4.5, 3.5, 3.0 


26.5 


5 


3 


4 


6.0, 5.0, 6.0 


17.0 



32 



3.3. timation of Par^eters 

t I 

It is convenitent Co use a and g from the pilot 
sample as initial guesses in the iterative technique for 
solving Eqs. 4 and 5 for the final estimates of a and g. 

The results obtained were; 

I * ’2 2 

a - 0,591 minutes, g » 0,590, 6 = 0.599 minutes , 

t 

6 = 0.774 minutes, 

= 0.591 minutaS) W2 = 1.066 minutes, 

= 1,922 minutes, = 3.468 minutes, 

w^ = 6.256 minutes. 

The experimenter is now in a position to test assumption c 

in section 3,1, that the deviations of yi^ at each value of 

o 

have a common variances 6 , For this the author will use 
Bartlett's test of Homogeneity of Variance, 

For this test. 



where 

K == 2.3026 [( df^)log - (Edf^ log 




1 

3(k.l) , 



k = 



and 



the number of populations of w's 
2 ssj^ + SS 2 + ... + ssj^ 

SdE 

^ X • 






In this example. 



L 



— — () Q3 

3(5-1) • ’ 
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= 8.97/30 = 0.299, 
P 

log = -0.524, 

P 

Zdf- log = —16.15, 



K = 2.3026 (-15.72 + 16.15) = 0.99, and 
y} = .99/1.08 = 0.917. 

Sincere? 95 4 ~ 9.48773 one cannot reject HO: ~ ^2 ~ 

2 * ’ 

••• = 6 ^ and the assumption seems to be valid. 

It is also possible to calculate 957* confidence 
limits on 



f 

E(u.) = w. * (t 55 .l>-^ "I i dj,. 

’ 1 



For the Circle K data, 

dj^ = 0.52, = 0.52, 

do=0.48 d=0.56, and 

4 

d_ = 1.31. 

5 

Since the experimenter is using the same estimate of 
experimental error for each array, the length of the confi- 
dence interval depends only upon n^^. The graph of vs 
is shown in Fig. 10 along with the 95% confidence at each w^. 
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Number of Vehicles in Food Herkot 
FIGURE 10. 

Time in Food Market (w) vs Number of Vehicles 
in Food Market (w) 
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3.4, LValuation of the riodel 

The author has estimated a, 3 and 6 for the traffic 
sink. It is now desired to examine the goodness of fit of 
the model to the physical system, j;xperimental error (6) is 

t 

always an indication of goodness of fit. The value 6 =C,774 
minutes may be good, it is impossible to say how good how- 
ever , since there is nothing to compare this with. As X'jas 
stated earlier, what is really hypothesized when it is 
assumed that c:(wj^) = ae^i is that the rate of change of w 
(as well as w Itself) increases as N increases. This assump- 
tion can be tested by comparing the model with a linear 
regression model. 

The linear model is C(wj^j) = a + 3Nj^ + ^<^ij 

where is the observation from the i^^ population, 

is the deviation from linearity of the i^^ 

population, and c is tho experimental error of the i 

population. Again it is assumed that the experimental error 

at each value of is normally distributed with the common 
o 

variance 5". This model has the advantage of testing 
directly the hypothesis that the deviations from linearity 
are equal to zero,^^^^ The analysis of variance for this 
model, using the same data is given in Table 3, 

Since E (MS departure from linearity) = 6^ + (k-2)*^+^Z^ ^’i^i 
and E(MS Residual) - 6^ we can test HO:Dj^*D 2 =* * * *D|^=0 v;ith 
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TABLE 3. Analysis of Variance for Linear Model 



Source 


df 


SS 


M.S 


Linear 

Regression 


1 


2^^ - 9 

b S m. (N,-N)'^=52.26 
i=l ^ ^ 


52.26 


Departure froTi 
Linearity 


5-2 = 3 


k _ = _ o 

S n.(w -w-b(N -N) =5.16 
1=1 1 i i 


1.72 


Residual 


35-5=30 


Z (w. .-w = -?.97 

1=1 i=l ^ 


0.299 


Total 


34 


k n 4 BE 9 

E (w. .-w) = 66.39 

i=l j = l i-J 


55.54 



the ratio 

7 ★ IS Departure from linearity 
MS Residual 

The results are N = 1.69, w = 2.56, 6 = 0.92, B = 0.97, 

Wj^ = 1.25, W 2 * 1.87, =* 2.33, = 3,79 and = 5.67, 

F = llg ~ 5.75, and ^ ^ 30 ^ l7l. Hence HO is 

rejected, and the author concludes that the linear model is 
not adequate for predicting E(w^), It has not been shown 
that S(wj^) = ae^^i is the best possible model for a traffic 
sink of the type considered. But it can be concluded that 
it is better than the linear model. The linear regression 
line is shown in Fig, 10. 
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If the exponential model had not given a better 

prediction of E(wj^) than the linear model, the experimenter 

would be forced to change the model of the traffic sink. 

The author feels the polynomial regression would do 

well in this situation since the experimenter has the 

ability to step the degree of pol3moraial until he finds 

the best fit of the .polynomial models. 

At this point, it can bo concluded that the expected 

waiting time of vehicles entering the traffic sink in 

6N 

question can be represented adequately by E(wj^) = ae^ 




I 

4 

4 
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CHAPTER 4 



COKTROLLU'JG TRAFFIC SHsTC CONDITIONS 



4.1, Introductio n 

After estimating a and 3 and satisfying himself that 
he can predict adequately waiting times for a particular 
8 ini:, the experLmenter still must answer one more question 
before the philosophy of Chapter 2 can be applied to a 
tunnel-sink system. Let Yj represent the expected number 
of arrivals/minute at sink^ of vehicles which have traveled 
through the tunnel. This value y., can be estimated from 
the origin and destination data available for the tunnel and 
must be the average over a certain time interval within 
(t ^ Again letting X. equal the expected number 

of arrivals per minute from all sources other than the tunnel 
for sinl:. the question is then, does the tunnel supply a 
large enough proportion Yj/(Yj+Xj), of the expected sink 
traffic to ensure that traffic controls at the tunnel will 
be reflected in the sink? 

If 907o of the traffic in the sink enters via the 
tunnel, then the experliaenter knows he has sufficient con- 
trol over the sink to employ the method of Chapter I, How- 
ever, if only 17o of the sink traffic comes through the 
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tunnel the eicperlnenter can do almost anything to the tunnel 
without affecting conditions in the sink. In other words, 
there is a critical value k* of k above which the experi- 
menter does not have sufficient control of the sink condi- 
tions to merit controlling the txinnel traffic. For k>k* ^ 
it is expected that the optimum waiting time in New Jersey 
will always be zero. This critical value of k must be 
estimated separately for every sink and depends upon a, 3, 
the expected operating conditions of the sink, and the 
expected rush hour conditions in the tunnel, 

Ilstim^tion of the Critical Value of k 

In order to estimate X*, it will be assumed that 
f(HQ,t,tj^), the function v^lch predicts the arrival of 
vehicles which have left the tunnel but have not yet arriv- 
ed at the sink, and g(N,t,tj^), the function which predicts 
the dissipation of vehicles initially in the sinl^, are 
knovm. Again let be the nximber of vehicles vdilch have 
entered the tunnel but have not yet arrived at the sink and 
let be the nimber of vehicles presently in the sink. 

A reasonable first approximation for f and g is to assume 
that they decrease linearly from their initial values. If 
H is depleted linearly and reaches zero in minutes and 
Nq decreases linearly to zero in minutes then; 
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where k 



1 




and 




v;here 0 = exp 1 acted flow for the next U 2 min, , 

fT = expected no. of vehicles in the tunnel for next 
min. , and 



Consider the second equation. The term (M^-rl^tj^/U^) simply 
estimates the nxambor of vehicles bound for the sink x»>hich 
have not yet arrived. If tj^ = 0, this term is and if tj^= 
U 2 this term is zero. Since t^ is assumed to be a constant. 
Vehicles enter the traffic sink with the same flox^ with 
which they exited the tunnel. If all vehicles were orig- 
inally in the tunnel, then (K -*i t./U-)/M indicates the per- 
centage of vehicles in the tunnel with this sink as their 
destination. Since 0 veliicles/mlnute exit the tunnel, then 
(ilQ-HQtj^/U 2 ) (^/FI) vehicles exit the tunnel each minute for 
this particular traffic sink, and t^ minutes later they 
enter the sink. Then, in general. 




w = ae 
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The term Is introduced to account for the time delay 
between exiting the tunnel and csntorins the traffic sink. 
Notice that the vehicles uhich enter from sources other than 
the tunnel are never removed from the sink. This implies 
that, besides the first-in, first out assxmption, it must be 
assumed that these vehicles will not depart from the sink 
before tj^ reaches its maximum values. Although these 
assximptions seem reasonable, they are not strictly necessary 
as will be shown later, 

4.3, 4^ j pple of Estimation Procedure s 

To deraonstrate the method of determining X*, it is 
hypothesized that a scaled down version of the Lincoln 
'funnel supplies the Circle K with veliicles. The Circle K is 

O 

used as a, p and 5 have already been estimated for it. It 
is assumed that t 3 = 1 rainute and the upper limit on tj^ is 
five minutes, i’ig. 11 shoxj^s the plot of E(w) vs tj^ for 
various values of X. Time in the tunnel (^ 2 ^ simply a 
scaled dox^m version of Fig. 6. x^ithout shock waves. The 
values of tj^ x^hich minimize T are shoxn in Table 4, For 
X > 1,4, the optimal waiting time (tj^) is always zero. This 
is because for X>1,4 the experimenter cannot exert suffi- 
cient control over the traffic sinJc to make waiting xTorth- 
xrhile. Uncontrollable soxirccs overshadow all the 
experimenters efforts to reduce the number in the traffic 



Time in Traffic Sink (Minutes) 
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FIGURK 11 

Time in Traffic Sink (w) vs Waiting Time in 

New Jersey (tj^) 

sink. Again using Fig, 6 as a guide, a similar analysis 
was made assuming that a shock vrave occurred in the tunnel. 
The results were the same and indicated X* - 1.4, 

In order to circumvent the first- in-first-out assump- 
tion and the assumption that no vehicle which enters from a 
source other than the tunnel vrill exit the sini< before 
tj^ =! the method of depleting N must be modified. It 

can be assumed that a constant percentage of the vehicles 
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TABLE 4, Linear Depletion of f and g 



X 


t]^ optL^n 


T mlnimura (minutes) 


0 


3 


8.9 


.1 


3 


9.4 


9 

• 


3 


9.9 


.3 


3 


10.6 


.4 


4 


11.2 


.5 


4 


11.9 


.6 


4 


12.9 


.7 


4 


14.1 


.8 


5 


15.4 


.9 


5 


17.1 


l.O 


5 


19.4 


1.1 


5 


22.5 


1.2 


5 


26.6 


1.3 


2 


30.3 


1.4 


0 


31.0 


1.5 


0 


31.0 



utiich are in the sinlc at the beginning of one minute will 
have e:-:ited by the end of that minute. At the beginning of 
each nijiuto, vehicles enter the traffic sinlc from uncon- 
trollable sources and possibly from the tunnel and combine 
with Vehicles already in the sink. By the end of this 
minute a certain percent of the vehicles have departed. 
Besides the previously mentioned advantages, this approach 
does not stipulate the manner in which vehicles leave dur- 
ing a given minute. This method does require an updating of 
N at the end of each minute before vehicles enter for the 
next minute, but still falls under the category of assuming 
the form of f and g. 
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?ig 12 shows v; as a function of tj^ if no shodc waves 
are encountered in the tunnel and if the constant percent 




Waiting Time in New Jersey 
(Minutes) 

FIGURE 12. 



Time in Traffic Sink (w) vs Waiting Time in 
New Jersey (t^^) 

depletion method is used to update N, For 1 > 0,6 it is 
seen (Table 5) that the experimenter does not have suffi- 
cient control over the sink traffic to make waiting before 
entering the tunnel profitable. Similar results wore 
observed with a shock wave in the txannel. T minimum was 
increased but X* remained 0.6, 

Ihe experimcjntor must now decide, from his know- 
ledge of the traffic sink, ^vhich method to employ in 
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TABLE 5. Constant Percentage Depletion 
of f and g 



X 


tj^ Optimum 


I minimum (minutes) 


0 


3 


8,3 


.1 


3 


8,9 


.2 


4 


9,8 


.3 


4 


10,6 


.4 


4 


11,8 


.5 


5 


13,7 


.6 


0 


15,9 


.7 


0 


15,9 



estimating X*. For example, the author feels that the first 
in, first out assumption is not unreasonable for the Circle 
K sampled. However, he does not believe that the first 
method can be realistically applied. While making a deci- 
sion for one vehicle, some of the contontrollable vehicles 
will enter and exit while 0 < tj^ < tj^ Although this 

will not be true if w is expected to be larger than , 

the author feels that for the Circle K, X* = 0,6 is a more 
reasonable limit than X* = 1.4, 

It is not difficult to estimate X, I/hlle sampling to 
estimate a and g , the easiest method is to record the 
initial value of K and then mark the time of each arrival 

t 

and departure. For this X can be determined. Precision 
requirements on X may necessitate drawing a larger sample 
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than is indicated in Chapter 2, however no theoretical 
problems are involved. 

If I < A.’', the experimenter can feel that he has 
enough control over the sink traffic to be able to reduce 
T by controlling tj^. 



CHAPTER 5 



SUi-itAKY MID CONCLUSIONS 



5,1. Sumniary 

The hypothesized nodel for predicting total time in 
the system is T = tj^ t 2 + constant + ae^^C^l). As Tj^ 
increases, N and hence E(w) decrease. The objective is 
to find the value of tj^ for which the reduction in T 
attributal to descreases in E(w) is offset by the increases 
in T caused by the growth of t^^ and 

Before attempting to use this model on an actual 
tunnel-sinlc systsn, two basic questions must be answered, 

1, Can the total time in the traffic sink be 
predicted adequately? To answer this question, the experi- 
menter must sample the sinks to estimate the parameters a 
and B, using, for example, proceedures outlined in Chapter 3, 

2. Does the tunnel, over which the experimenter has 
control, supply a large enough percentage of the input to 
the sink so that the reduction of N will overcome the 
increases in T caused by tj^ and tp7 This problem is dis- 
cussed in Chapter 4 and involves estimating X. and X*. 

For the tunnel, it was fovind that t 2 increases with 
tj^ as in Fig, 6 of Chapter 2, For the Circle K it was found 
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that a = 0.59, t! = 0.591, \ =0.47, and X* = 0.6. Thus, if 
a scaled down version of the tunnel supplied vehicles Co 
the Circle K, the optimum would be 4 minutes and would 
result in an expected total time in the system of 13,3 
minutes. These results are summarized in Fig, 13, 



T 




FIGURE 13. 

Total Time in System (T) vs Waiting Time in 
New Jersey (t^) 

Since it is reasonable to believe that t^ will actually 
decrease as tj^ increases and thus help decrease T, the 
assumption t^ = constant is a conservative one. 
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5.2. Acplication of the i ^ odel to the Lincoln IXin nel System 

There are several stxirabllng blocks which must be 
overcome before the idea can be applied successfully to the 
Lincoln Tunnel - New York City system. 

1. There are no vehicle waiting areas available in 
New Jersey, 

2. As of this writing, the authorities had not yet 
set up a traffic control in the tunnel so as to avoid shock 
waves . 

3. Because of the sizes of the traffic sinks, the 
number of vehicles involved, and the large number of sink 
entrances and exits, the sampling problems in New York are 
expected to be very difficult, much more so than for the 
example developed in Qiapter 3, 

Due to the myriad of entrances into any New York 
sink the question of controllability posed in Chapter 4 is 
very real. It is unlikely that the Lincoln Tunnel supplies 
a large enough percentage of the traffic into the Central 
Business district of New York City to make this method of 
control applicable to all 24 traffic sinks. 

5. 3. Extension to U. S. Naval Ships 

The theory and methods discussed in Chapter 2 apply 
to any transportation systan in which controlling the 
vehicles as they enter will affect the conditions within 
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the system that these vehicles will encounter. To the 
author, one of the most obvious extensions of these control 
procedures Is to Naval ships entering port thirough the con» 
stralnts of a channel, inany of the problems encountered In 
New York are alleviated In this channel-harbor system and 
some new difficulties arise, but the basic control theories 
remain the same. 

Tlio traffic sinks within the harbor are the deep 
draft piers, the shallow draft piers, fueling piers, 
anchorages, possibly ammunition loading areas and possibly 
several points within a shipyard. There are fewer sinks 
than in New York and they are more distinct and easier to 
define. The sampling procedures are simpler than those in 
New York and probably wovtld require little more than sampl- 
ing the watch officers' logs of selected ships. 

There would be no requirement to construct waiting 
areas at the channel entrance, and the problem of control- 
lability Is completely removed. The Senior Officer Present 
Afloat (SOPA) has absolute control over all vehicles which 
would enter the traffic sinks. Thus \ would be zero and X.* 
non-existent. Because of the speeds and spaclngs Involved 
it is much easier to predict conditions In a channel than 
to predict conditions in a tunnel. 

The interactions between ships attempting to be 
serviced (to moor) In the same sink are just as real and as 
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bothersome as the delays caused by trucks trying to off load 
cargo in the same area. It is believed that because of the 
simplifications inherent in a channel-harbor system, the 
control procedures outlined in the preceding chapters can 
be applied successfully to these systems at this time, 

5, 4, Conclusions 

The model of Chapter 2 is by no means complete. For 
example, no mention has been made of the Holland Tunnel, 
another facility under the control of the Port of Now York 
Authority, 

The initial study of a large problem, must by 
necessity, involve many simplifying assumptions. These 
enable the experimenter to isolate major problem areas as 
was done in Chapters 3 and 4, 

In conclusion, it is believed that the idea of 
reducing traffic delays in New York City by controlling 
vehicles before they enter the tunnel is possible and 
worthwhile. Authority officials feel that it will be years 
before they actually begin experimenting with the ideas 
developed in this thesis. But the problem of traffic con- 
gestion , outlined in Chapters I and 2, must someday be 
faced. It is hoped that this study constitutes a first 
step. 
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D 1 ‘^FNS ION TCDNJ ( 2 ) ^ T NW TO ( 2 ) , 5> T I T ( 2 ) , AEND ( ? ) « TQDSP ( ? ■» 

1 100 ) . T IS ( 2 ) ( 2 . 1 00 ) , ^STOr ( 2 I , T-^AOS ( 2 , 1 '"0 ) , OB ( ? , ) 

2,T 10(2) »TT( 2) .TODS7( 1'‘'0) .TOaJ7(10C) ,BM(?5»5) » 

IP'^ ( 25 »S ) *S : ( 100 ) .07 ( lOG ) , ENTER ( 2 ) . 1 COL ( 25 ) » -V | ♦ 
4AA(2 ) .&( 2 ) 

N'ROV = 0 
\rOL = 0 

READ eOO » HOUR . TN,RHO ,RTTO 
SOO FORMAT ( I 1^, AF10.5 ) 

PFAD pot. (TppnJ(JJ), JJ= 1 ,'M 
pDl FOPk/AT ( PFI '■',5 ) 

RFAD , flNWTG(JJ), JJ = I 

READ aO] » ( ^ T I T ( JJ ) ♦ JJ= l,v) 

REAI'' PO] 1 ( aEmD( JJ ) , JJ = l,v) 

Rtad . ( aSTOP ( JJ 5 » JJ = 1,M) 

REAP 8C1»(TT!JJ), JJ = 1,M) 

READ 801- (ENI'erijJ) , JJ - 1,V) 

DO 802 JJ = 1 
I A = AEND(JJ) 

602 read 801 . ( TuDSR ( JJ* JJJ) .JJJ 
DO JJ = 1 

= AFVP(JJ) 

pO-a RFAD POT * ( SB ( JJ * JJJ ) ♦ JJJ = 

DO 8'" A JJ = 1*M 
I A = A<^'T^'^(JJ' 

P'^h read poT , ( TOAGR ( JJ * JJJ ) » JJJ 
DO 805 JJ = 1*M 

1 A r.- ASTO^^' IJ) 

805 Reap 901 * ( jo ’--jJ > * -*oj = 

READ 801 ♦ ( V.'G! I ) * I = 1 *M> 

READ 8 1 ♦ ( A A ( 1 ) , I = 1 , M ) 

'^^AD 80 1 * i P ( I ) « 1 = 1 *M ) 

2 = '^.o 

DO 890 I - 1 V 25 
A ( 1) - 7 
1 7 + 1.0 

TIT = 0 . 

O'"* 805 JJ = 1 *M 
p06 TIT = TIT + STiTvJJ) 

DO 007 JJ = 1 *M 

I SINS = AEND(JJ) 

r.o P07 JJJ = ItlSPO^ 

8O7 TIS(Jj) = TIS(JJ) + SB(JJ*JJJ) 

DO 808 JJ = 1 *5 
I'Mnk = A'‘T0P(JJ) 

DO POP JJJ = 1 * ISINK 
p^'p TtO(JJ) = TIO(JJ) 4 0o(JJ*JJJ! 

I A = 0 

3 38 IF ( T N - 8 0 . 0 ) 5 1 0 *''11*511 
511 HO'JR = HOUR 4 1.0 
TN = 1.0 
60 TO 51A 



= i ♦ I A ) 
1 . I A ) 

- ] . I A ) 
1 . lA) 
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5 1C IN = TN 4 2.0 

RT^j = HO'JR 4 T .\/6 0,0 

ir (PTN - PTTQ) 994P ,9949 .'^^•'-•8 

9040 poir^y , hOOR , IN' 

995 c FCR'^ATf// 48X,15 HThE TIVE 10 NOW, 2 X , F? . 0 , 3X , F3 . 0 ) 
NT IM.O = 0 

DO 19 I = 1 
I OINK = AFND( I ) 

DO lO J = 1 , IS INK 

7 K = J 

4 IF (T0DSB(I*<1 - 0.0) 2.2*3 
? IF (K - I 9 INK) 9,11,11 
9 < = K + 1 

GO TO 4 

3 IF (J - K)9,S,9 
9 Torso ( I , J ) = T0DS8( I ,<) 

T''DSBt I .K ) = 

F ir (TOP^P(T,J) - PTN) 6,6,1'^ 

0 T IS( I ) = TiS( ! ) - S 8 ( I ,K) 

AFND ( I ) = Af Nir ( I) - 1.0 
SfUl.K) -4 G.O 
TCDG3M.J) = O.C 
GO TO 7 

10 CONTINUE 

11 I INK = AFND ( I ) 

DO 19 J = 1 , I SI NK 

K = J 

13 IF (Snil.K) - (.0) 12,12,17 

12 v' = V" + 1 

GO TO 13 

17 IF (J - K) 24,18,14 

14 FP ( I , J ) = SB ( 1 ,K ) 

SR M ,F ) = 0.0 

19 CONTINUE 

DC 40 I = 1 ,M 

rrwp = o.F 
ICI^.K = afTO?<I) 

DO 30 J = 1 , I S I NX 

27 K = J 

24 IF (ICAOBfl.K) - O.J) 22.22,23 

2 7 Jr(^' - T'^P’K) 1,31 

2 F K = K + 1 

TO 24 

23 IF {J - K) 29,28,29 

2P TOAOBII.J) * rOAGR(I,<) 

TOAQB ( ! ,K ) = O.C 

28 IF (TOAOF(I.U) " RTN) 26,26,30 

26 T I 0( I ) = T ID ( I ) - QB ( I »K ) 

; /^TOP ( I ) = Ac roP( I ) - 1.0 

T fwp = ‘-’P 4 OB ( I ♦ K ) 

CR(!,<) = 0.0 

T^AOR ( I , J ) ~ 0.0 
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G T O 2 7 




30 


C CMT I NOT 




31 


ISINF = AST 


OP ( I ) 




IF fISINS - 


0) 39,39,100 


O'"' 


on TO J = ) 


. I SI NS 




F r J 




3 3 


IF (OB ! I ) 


- 0.0) 32.32.37 


3? 


K = < +1 






GO TO 33 




37 


IF (J - F ) 


34 ,38.34 


34 


0‘3( I ,J) = OP. ( I ,K1 




03(1.'^) rs 0 


.0 


39 


0 r>N’ ’■ I on ifT 
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J ^ ArNDI I ) 


^ 1 .0 



S^M . ' > = TilMP + ENT PR ( I ) 

T 1 S ( 1 ) = M S ( I ) + SR ( I « J ) 

•V T = A A ( 1 ; « e X F ( o 1 I ) * T I s { 1 1 ) 

TOr/PR(I,J! ^ RTN + WT/6C.0 
AFNr>(I) = J 
AO CC'NTINOT 

DO ‘=>0 1 -■ 1»M 
J = A STOP ( 1 ) 1- 1.0 
QR(I.J) = ( ST I T 1 I) / T I T ) #PH0 

T'^AOB ?NJ) = RTN + TT(I) 

I IO( I ) = T I 0( I ) + on( I , J ) 

AS TOP ( I) = J 
T i T = TIT - OB ( ! . J ) 

SO ST I T ( I ) i S T I T ( I ) - GFM I , 

512 RHO = -01.6 + 2C.0 * RTN - ?.0 » RTN **2.0 
C^TETR = RnO 
AGTV = 0.0 
DO AOO I ^ 1,M 

IF (TODNJ(I) - 0.0 ) AOO.AOo.AGl 
AOl ir (TOPNJ(I) - RTN) AO? ,402. 

40? STIT(I) = “^TITll) + TNWTG( 1 > 

TI-^ = TIT + TNWTG(I) 

AGTf/ = aoTM -t- TNV.'TG(I) 

PC' IN'" 9951 . TNWTG( I ) , I 

9951 FORMATf 1X.F5.0,35HVEHICL£S OEPARTE'D N.J. FOR S I .NF N'.JO 
14HBER, 15; 

TODNJ (I) = 0.0 

TNWTG(I) = O.C 
400 CONTINUF 

IF (AGTT* -• DTlTM) 419,333,333 
419 DTFTN = DTFTM - AO-TM 
AGT^^ = "'.O 

4 24 NT IMS - NT IMS t 1 

ir (nUmS - V) 428.426.333 
42fl IF (lA - M) 422.421.421 

421 lA - 1 

GO TO 423 

422 lA = lA + 1 



A?? IF (TOONJ('A! - C.O) k2i**‘*cb 
GO TO 70^ 

4 76 TOO''!J(IA' = RTM + (A(,NHOWl ) /60.0 
TNWTG (lAO ^ NCOL 

t^P I\T P957 »NCCL , A ( NRO'-: ' . 19^- ( 1 . 1 ) 

996? FORMAT ( 1 X , I 3 , ?^,HVhHia.rS WfPF TCLD TO V'AlT,c-6.0, 

11^F’‘INUTrs !N N . J , , 6X . PBHTCr IIMF 'JOTIL DFPAFnuRF V-.' 
?6MILL PF . I 6 *BHM!NUTF:s. ) 

IF (TCD'^J(IA) - RTM A?7,42 7,42A 
42 7 STIT(IA) = oTITflA) + TNwTC. (lA) 

TIT =T I T ■*■ T,\WTG( I A ) 

AGTV = AGTM -t- TNWTG(IA) 

PPPIT 995 M TNS'TG < I A ) , I A 
TOONJtlA) = 

TGWTC(IA) ^ O.C 
IF (AGTM - DTFT*^) 424,333t3?3 
7 or TIT7 = TIT 

P T 7 = 7 T V 

5 T I T 7 = S T I T ( ! A ) 

T2 = TTIIA! 60,0 
I END = ArvP(lA) 

I ATOP = A STOP ( 1 A ) 

DO 70 1 J = ]»!ENO 
S7(J) = FB(IA.J) 

7C ] T ODO ( J ) = TOPSB ( ! A , J ) 

DO 70? J = 1 » I SI CP 
Q7 ( J ) = OB ( f A » J ) 

70? Tr'AO-'(J) = TOAOB(1A»J) 

T I F 7 - I ' ( I A ) 

T 107 ^ T 10( lA * 

'^ = 0, 

DO 77 ^ II = ] .75 
DO 704 J = 1,5 

704 BM( II . J ) = S 

G = 5 + } . f' 

7C3 CONTINUE 

5 = 0.0 

DO 75 c II ^ 3,25 
T = PTN7 + S/60.0 

RHC7 = - 3i.6 + 20.0 *T - 2.0»T**2.0 

CON = ( .61 ) / 2 . » . C 043 ) + ( ( . 6 ] *» ( 2 ) - 4 .*, 0043 * ( . 44 - RH 07 ) ) 

1 /( 7 . », P 047 I 
?04 VFL = PHO '^ / CON 

T 1 = ( B 000, 0/52 39 , 0 ) /VFL 
DO 705 J = 1,5 

705 5M(II»J) = bM(Ii,J) + T1 + T2 
708 IT = 1 

200 RTN7 = PTN7 + 1.0/6(y.C 
I^INK r -. I END 
DO 210 J = 1 , i SINK 
?07 K - J 

204 IF(TODS7IK) - 0,0) 2O?,?02.?03 
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GO TO 

?03 IF(J - K) ,20a ,209 

20? TrD<^7(J; = T0DS7(<) 

TrDF7(K) = C'.n 

.’Oe IF (TCO'^TIJ) - RTn7) 206,206*210 
?<"6 T I ?7 = TIG' - S7 ( F ) 

I END = Il N'O - 1 
S7(<) = 0, 

T0D.S7(j) = 0.'' 

GO TO 207 
510 COMTINl'f^ 

2)1 lEINK = TEND 

pr\ .;>ip j = 1,1 5 INK 

K - J 

2)3 IF (f)7(K) - F.n) 2 1 2 , 2 1 2 , 2 '' 7 

212 K = K + 1 

GO TO 213 

217 IF (J - K) 219,218,2)9 
■■I 9 S7(J) = S7(K) 

S7 (K ) = 0.0 

218 COMTINUc 

irMP = 0.0 

I S INK = I STOP 
D<^ 230 J = 1 , ISINX 
7?7 ^ = J 

229 IF (T0A07{<) - 0.0) 222,222,223 

22 2 If- ( K - I S J NK ! 2 2 F , 2 5 1 , 2 3 1 

225 K = K +1 
GO TO 229 

22’ IF (J - K) 22 «, 2 2 3, 229 
229 T0A07IJ) = T0A07(K) 

TC907(K> = t-.O 

228 IF (T0A07(J) - RTN7) 226,229,230 

226 T 107 s T 107 - 07 ( K ) 

I STOP = IS TOP - 1 
TE"P = TEPr + 07(.<) 

07(0 - O.C 
T0A07(J) = 0.0 

GO TO 227 
23C CCNTINUG 

231 I SINK = I STOP 

If (I SINK - Z) 239,239,900 
90'* DO 238 J - 1,1 SINK 
K = J 

233 IF (07(K) ~ 0.0) 232,232,237 

232 K = K +1 
GO TO 2’3 

237 If (J - K) 239. 238, 2 ’9 
239 07(J) = 07 (K) 

07 (K! = 0,0 
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?38 CONTiNLir 
?39 J = rPND + 1 

:.7(J) = TF^'P ♦ C.\TFR(IA) 

11^7 = T ! ‘^7 + S7 ( j ) 

VT = AA ( ] A 1 :»r XPFF ( !') ( I A ) ■>■ T 1 ?.7 ) 

TOr'A7(J) = RT’i7 + WT/ 60.0 
If^■0 ^ J 
J = I T r. A' 1 

07(J) = ('-T!T(IA< / T;T7)*Pm07 

T0A07 ( J ) - KIN -t TT ( I A ) 

TI07 = T107 f C7(J) 
iSTOP = J 

ST I F7 = STITT - 07 ( J ) 

PH07 = -31.6 + ?;C.0T(R7,M7) - 2 . 0 » ( P T ’'M ^ . 0 

7A9 CONTINUE 
ADD = C'.O 
DO 7 P.0 J = 

TI57 = riS7 f ADD 

WT = AA ( I A ) *rxPFF ( n ( : A ) it r I 7 ) 

3‘' ( n ♦ J ) = BP ( I 1 . J ’ + WT 

72^ AOi^ - ADD + 1.0 

5 ^ S 4 1.0 

7 50 T IS7 = T I :,7 - ADD 
DO 950 II = ) *75 
DO 950 J = 1*5 

95C IRP(1I*J) - BV(1I,J) + 0.5 
DO 9950 I = 1*75 
9961 FOPMAT (1 X , 5 I 6 5 

9960 PRINT 996], (IE'P(],J). j = T, , 5 ) 

DO 962 II = 1*7^ 

I rOL (II! - ] 

DO 952 J = 7*5 

IF - I Of-' : I 1,1)) 06-^ .OCT ,957 

qr.T I PM ( M , 1 ) c TPM ( T I , J ) 

I COL (II) = J 
95 7 CONTIK’L"=' 

NCOL --- ICOL(l) 

NPOW = 1 

DO oc/, T 1 = 7*75 

IF ( I 9 V ( II , 1 ) - I P ■*! 1 1 , 1 ) ) 9 5 5 , 9 P 5 , 9 5 c. 

95 5 1 f P ( 1 ♦ 1 ) = IBP ( I I * ] ) 

NOOL = I COL (II) 

MPO'’' - I I 
95a con TIN"" 

C,o rc A? 6 

99AP STOP 
FND 
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